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PREFACE 


The FORT RAN *-63 language contains all of the features of its predecessor, 
FORTRAN-62 and forms an overset of the FORTRAN II language. The FORTRAN-63 
compiler adapts current compiler techniques to the particular capabilities of the 
Control Data ** 1604 and 3600 computer systems. Emphasis has been placed on pro¬ 
ducing highly efficient object programs while maintaining the efficiency of compi¬ 
lation of FORTRAN-62. 

This reference manual was written as a text,for FORTRAN-63 classes and as a 
reference manual for programmers using tfie FORTRAN-63 system. The manual 
assumes a basic knowledge of the FORTRAN language. 

Reference material consists of three volumes: 


Volume I Arithmetic and Logical Information 


Volume II Input - Output 

Volume in Will contain programming instructions for type 

non-standard arithmetic and instructions for 
compilation and execution of FORTRAN-63. 


*FORTRAN is an abbreviation for FORmula TRANslation and was originally developed for 
International Business Machine equipment. 


Ill 


**Registered trademark of Control Data Corporation 
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ELEMENTS OF FORTRAN-63 


1.1 

CHARACTERS 


1.2 

OPERATORS 


This chapter presents the character set used by FORTRAN-63 and demonstrates 
how entities within it are used to form the identifier. Arithmetic, relational and 
logical operators are listed along with their meanings, and ideas relating to number 
are presented in terms of the fundamental definitions of quantity, variable and 
constant. 


FORTRAN-63 uses the following character"set; the conventional FORTRAN 
definitions apply: 

The alphabet A through Z 

The Arabic numerals 0 through 9 

The special characters + - = / ( ) . , $ * space 

The special character $ is a statement separator which may be used to write more 
than one statement to a line. 


The operation symbols used in replacement and conditional statements are 
tabulated below: 


Mathematical 

Symbol 

Meaning 

FORTRAN-63 

Operators 

Classification 

+ 

Addition 

+ \ 


- 

Subtraction 



4- 

Division 

/ 

Arithmetic 

X 

Multiplication 

* ' 


[ ]" 

E xponentiation 

** 


= 

Equal to 

■ EQ. ^ 



Not equal to 

.NE. J 


> 

Greater than 

.GT. ( 

/ 

i Relational 

> 

Greater than or equal to 

.GE. 1 


< 

Less than 

.LT. ' 


< 

Less than or equal to 

.LE. j 


A 

Conjunction 

.AND. 

> 

V 

Disjunction 

.OR. j 

> Logical 

- 

Negation 

.NOT. 

I 
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1.3 

IDENTIFIERS 


1.4 

QUANTITIES AND 
WORD STRUCTURE 


Mathematical 

Symbol 

Meaning 

FORTRAN-63 

Operators 

Classification 

n 

Logical product 

.AND. \ 


U 

Logical sum 

.OR. > 

Masking 

~ 

Complement 

.NOT. / 


<- 

Is replaced by 

- 

Replacement 


Identifiers in FORTRAN-63 fall into two classes: numeric and alphanumeric. The 
numeric identifiers are: 

Bank designators in the 3600 system represented by a single octal digit. 

Statement numbers represented by a number appearing in card columns or 
printer positions 1 through 5. This number is in the range l^N—99999. 

Block COMMON identifiers. 

The alphanumeric identifiers which name variables, arrays, subroutines, functions 
and the like may be from one to eight alphanumeric characters; the first of which 
must be alphabetic. Spaces in any identifier are squeezed out; A aa 6 is the same 
as A6. 


Examples 


A156 

SINEX 

ALPHA 

HEGEMONY 

G 

LUX31Z 

HERA 

A1B2C3D4- 


FORTRAN-63 manipulates floating point or integer quantities. Floating point 
quantities have an exponent and a fractional part. The following classes of numbers 
are floating point quantities. 


REAL Exponent and sign 11 bits; fraction and sign 37 bits; range of number 

(in magnitude) 10~ 308 — N — 10 308 and zero; precision approximately 
11 decimal digits. 


DOUBLE Exponent and sign 11 bits; fraction and sign 85 bits; range of number 
(in magnitude) HP 308 — N — 10 308 and zero; precision approximately 
25 decimal digits. 


COMPLEX Two reals as defined above. 


The following classes of numbers are integer quantities: 


INTEGER Represented by 48 bits, first bit is the sign; range of number (in 
magnitude) 0 — N ^ 2 47 -1; precision is up to 15 decimal digits. 
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1.5 

CONSTANTS 


1.5.1 

INTEGER CONSTANTS 


1.5.2 

OCTAL CONSTANTS 


1.5.3 

HOLLERITH CONSTANTS 


LOGICAL 1 in bit position 47 represents the value TRUE. 

0 in bit position 47 represents the value FALSE. 

HOLLERITH Binary coded decimal (BCD) representation treated as an integer 
number. 

A FORTRAN-63 program may contain any or all of these classes of numbers in the 
forms of constants, variables, elements of arrays, evaluated functions and so forth. 
Variables, arrays and functions are associated with types assigned by the 
programmer. The type of a constant is determined by its form. 


To define constants let: 

n be a string of decimal digits 

s be a scalar with a maximum of three decimal digits 
o be a string of octal digits 
h be the length of a Hollerith field 
f be a Hollerith field 
R be a Real 


n denotes an integer whose range, precision, et cetera, is as defined above. 

Examples 

63 3647631 

247 2 

314159265 464646464 


Octal constants may consist of up to 16 octal digits. The form of this constant is oB. 


Examples 

7777777700000000B 
7777700077777B 
7777777777777700B 


2323232323232323B 

77B 

same as .NOT. 77B (section 2.7) 
or -77B 


A Hollerith constant is a string of alphanumeric characters of variable length of the 
form hHf, where h is an unsigned decimal integer between 1 and 8 representing the 
length of the field f. Spaces are significant in the field f. When h is less than 8, 
the representation in the computer word is left-justified with BCD spaces filling 
the remainder of the word. An alternate form is hRf. When h is less than 8, the 
internal representation is right-justified with zero fill. 
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Examples 


1.5.4 

FLOATING POINT 
CONSTANTS 


6HC0GIT0 

4HERG0 

3HSUM 


8RCDC 3600 
8R ** 
1H) 


REAL Real numbers are represented by a string of up to ten digits. A real 

constant may be expressed with a decimal point or with a fraction and 
an exponent representing a power of ten. The forms of real constants 
are: 

+ 4 - + + 

n, n n, . n nE- s . nE - s n. nE - s n. E - s 

The plus sign may be omitted for positive s. The range of s is 0 ^ s 
^ 308. 


Examples 


3.1415768 

31.415 92 E-001 

314. 

.31415 E01 

.6749162 

.31415E001 

3 14159E-05 

.31415E+01 


DOUBLE Double precision constants are represented by a string of up to 

25 digits. A double precision constant has forms analagous to the 
forms of reals, with the E replaced by D. The forms are: 

4 " 4 " 4 " 

n.nD n.D . nD nD - s n. nD - s n. D - s 

The plus sign may be omitted for positive s; the range of s is 
0 ^ s — 308. The D designator must always appear. 


Examples 

3 .1415926535897932384626D 31415.D-004 

3.1415 D 379867524430111D+001 

3.1415 DO 
3141.598D-03 

COMPLEX Complex constants are represented by a pair of reals enclosed in 
parentheses with the reals separated by a comma: (R 19 R 2 ), R 1 
represents the real part of the complex number and R 2 represents the 
imaginary part. 


Examples 

The complex numbers 1. + 6.55i, -14.09 + .0001654i, 15. + 16.7i, and 
-i are represented in FORTRAN-63 as: 

(1., 6.55) (-14.09, 1.654E-004) 

(15., 16.7) (0., -1.) 
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1.5.5 

WORD STRUCTURE 


The word structure of the quantities in FORTRAN-63 is shown below: 











1.6 

VARIABLES FORTRAN-63 recognizes two kinds of variables, each is represented by an 

alphanumeric identifier. A simple variable represents a single quantity; a sub¬ 
scripted variable represents a single quantity within an array of quantities. The 
identifier appears with a subscript list enclosed in parentheses. The subscript 
list has the form (S 1 ) or (S 19 S 2 ) or (S 19 S 2 , S 3 ) where Sj may be standard or non¬ 
standard. 


1.6.1 

ARRAY STRUCTURE 

AND SUBSCRIPT FORMS Elements of arrays are stored columnwise in ascending order of storage location. 

A 3 by 3 by 3 matrix illustrates the storing process: 



The planes are stored in order, starting with the first, as follows: 


—*-L 

a —*-L+3 

a , 1 -L+24 

111 

121 

133 

—-L+l 

& „ —L+4 . 

.. a —*- L+2 5 

211 

221 

233 

— L+2 

a 00 „—»-L+5 

a —*- L+2 6 

311 

321 

333 


The maximum permissible number of subscripts appearing with a variable is 
three. The structure of the subscript is flexible within the classes standard and 
non-standard. 


1.6.2 

SUBSCRIPT RULES SUB1 A standard subscript has one of the following forms: 

C 

C*m 
m - d 
m 

C*m - d 
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where C, d are unsigned integer constants and m is a simple integer variable. 

SUB2 A non-standard subscript is any arithmetic expression used as a subscript, or a 
subscripted subscript. 

SUB3 The location of an array element with respect to the first element of the array is 
a function of the array dimension, type (3.1), and the subscripts appearing with the 
array identifier. In general, given DIMENSION A(L,M,N) the location of A (i,j,k) 
with respect to the first element A of the array is given by 

A + {i - 1 + L (j - 1 + M (k - 1) )} 

The quantity in braces is called the subscript expression. For standard subscripts, 
the subscript expression, when evaluated, is an integer. For non-standard sub¬ 
scripts that are arithmetic expressions, the subscript expression is truncated after 
evaluation. 


Examples 

1. Referring to the matrix in 1.6.1, the location of A (2,2,3) with respect to 
A (1,1,1) is 

Locn ( A(2,2,3) } = Locn {A(l,l,l) } + {2-l+3(l+3(2) )} 

= L + 22 


2 . 


Given DIMENSION Z 
The location, z, of Z 

z = Locn {Z(l,l,l)} 

= Locn {Z(l,l,l) } 

= Locn {Z (1,1,1)} 


(5,5,5) and I = 1, K = 2, X = 45°, A -7.29, B = 1.62. 
(I * K, TANF (x), A—B) with respect to Z (1,1,1) is: 

+ {2-1+5(1-1+5(4.67) )} Integer part 

+ {117.75} Integer part 

+ 117 


SUB4 FORTRAN-63 permits the following relaxation on the representation of subscripted 
variables: 

Given Array A declared with 3 dimensions as in DIMENSION A(D 1 ,D 2 ,D 3 ) 
where the Dj are integer constants. 


then A(I,J,K) implies A(I,J,K) 

A(I,J) 

implies A(I,J,1) 

A(I) 

implies A(I,1,1) 

A 

implies A(l,l,l) 

similarly, for A^^Dg): 

A(I,J) 

implies A(I,J) 

A(I) 

implies A(I,1) 

A 

implies A(l,l) 


The converse does not hold. The elements of a single-dimension array, A(D 1 ), 
may not be referred to as A(I,J,K) or A(I,J). 

Array allocation is discussed under DIMENSION, (3.2) and array structure (1.6.1). 
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1.7 

FORTRAN-63 

LANGUAGE 

STATEMENTS 


Examples 

Simple Subscripted Variable Subscripted Variable 


Variable 

(Standard) 

(Non Standard) 

FRAN 

A(I, J) 

A(MAXF(I,J,M) ) 

P 

B(I+2,J+3,2*K+1) 

B(J,SINF(J) ) 

Z14 

Q(14) 

C (I+K) 

ESTRUS 

P(KLIM,JLIM+5) 

MOT ZO (3 *K*ILIM+3.5) 

MAX3 

. SAM(J-6) 

WOW(I(J(K) ) ) 

I 

B(l,2,3) 

Q(l,-4,-2) 



The FORTRAN-63 elements in this section are combined to form statements, the 
basic program element. These statements are either executable or non-executable. 
An executable statement performs a calculation or directs control of the program; 
a non-executable statement provides the compiler with information regarding 
variable structure, array allocation, storage sharing requirements, and so forth. 
FORTRAN-63 statements are listed in Appendix C. 



§ 
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EXPRESSIONS AND 
REPLACEMENT STATEMENTS 


2 



There are three kinds of expressions in FORTRAN-63: Arithmetic and masking t 
expressions which have numerical values and logical expressions which have truth 
values. These expressions may appear in arithmetic, masking or logical replace¬ 
ment statements, in IF statements or in arithmetic function statements. 


2.1 

ARITHMETIC 

EXPRESSIONS 


All expressions are combinations of operators and operands. For the arithmetic 
expression, these entities are: 

Operators: ***/ + _ 

Operands: Constants 

Variables (Simple or subscripted) 

Functions (Chapter V) 

Two or more arithmetic expressions can be combined to form another arithmetic 
expression and so on. Parentheses are used to direct the order of evaluation of 
the expression. 


2.1.1 

RULES OF 
FORMATION 


A1 The hierarchy of arithmetic operations is: 


** 

exponentiation 

class 1 

/ 

division 

class 2 

* 

multiplication 


+ 

addition 

class 3 

- 

subtraction 



A2 Any variable (with or without subscripts) or constant, or function, is an arithmetic 
expression. These entities may be combined by using the arithmetic operators to 
form meaningful algebraic arithmetic expressions, subject to the following rules 
and definitions: 


t Called Boolean in earlier versions of FORTRAN. 
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1. Let op be an arithmetic operator and X,Y be arithmetic expressions. 
The form X op op Y is never legitimate. 

2. If X is an expression then (X), ( (X) ), et cetera, are expressions. 

3. If X,Y are expressions, then 

X + Y 
X - Y 
X/Y 

X*Y are expressions. 

4. Expressions of the form X**Y and X**(-Y) are legitimate. They are 
subject to the restrictions in Mode of Arithmetic Expressions (2.2.1). 


5. Implied multiplication is permitted, but only in the following four ways: 


constant (. . .) 
(. ..)(.. .) 

(. . ,)constant 
(. . .)variable 


implies constant * (. . .) 
implies 

implies (...)* constant 
implies (...)* variable 


A3 In an expression with no parentheses or within a pair of parentheses, in which 
unlike classes of operators appear, evaluation proceeds in the order stated in Al. 
In these expressions where operators of like classes appear, evaluation proceeds 
from left to right. 

A4 When an arithmetic expression contains a function, the function is evaluated first. 

A5 In parenthetical expressions within parenthetical expressions, evaluation begins 
with the innermost expression. 


Examples 

Expressions 

A 

3.141592 
B + 16.8946 
(A - B(I, J + K) ) 

G * C(J) + 4.1/ (Z(I+J,3*K) )*SINF(V) 

(Q + V(M,MAXF(A,B) ) * Y**2) / (G*H-F(K + 3) ) 

-C + D(I,J) * 13.627 

In the following examples R indicates an intermediate result in evaluation; it does 
not necessarily imply an object language store. 

A**B/C+D*E*F-G is evaluated: 


A**B - 



Rj/C - 

— R 2 


D*E - 

— r 3 


r 3 *f — 

— R 4 


R 4 +R 2 

— R 5 


r 5 -g — 

— R 6 

evaluation completed 
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A**B/(C+D)*(E*F-G) is evaluated: 

A**B -►Ri 

C+D -►- R 2 

e*f-g- R 3 

R 1 /R 2 -R 4 

r 4 *r 3 -*- r 5 evaluation completed 

If the expression contains a function, the function is evaluated first. 

H(13)+C(I, J+2)*(COSF(Z) )**2 is evaluated: 

Z - Ri 

COSF(R 1 )—R 2 

r 2 *r 2 -R 3 

R 3 *C(I, J+2)-*-R 4 

R 4 +H(13)—R 5 evaluation coinpleted 

The following is an example of an expression with embedded parentheses. 

A*(B+( (C/D)-E) ) is evaluated: 

C/D -R x 

R x -E -R 2 

r 2 +b- *- R 3 

R 3 *A -*- R^ evaluation completed 

A*(SINF(X)+1.)-Z/(C*(D-(E+F) ) ) is evaluated: 

SINF(X)-► R x 

Ri+l.- R2 

r 2 *a- r 3 

E+F -- R 4 

-R 4 *- r 4 

R 4 + d *- R 5 

r 5 *c- r 6 

-z - r 7 

R 7 /R 6 -+- r 8 

R 8 +R 3 -R g evaluation completed 

2 . 1.2 

INTEGER 

ARITHMETIC,CAUTION In both the 1604 and 3600 computer systems, dividing an integer quantity by an 

integer quantity always yields a truncated (least integer) result; thus ll/3 = 3. 
For this reason, plus the fact that expressions containing operators of the same 
class are evaluated from left to right, the expression I*J/K is not necessarily the 
same as J/K*I. For example, 4*3/2 = 6 but 3/2*4 = 4. 
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2.2 

MODE OF ARITHMETIC 

EXPRESSIONS FORTRAN-63 permits full mixed mode arithmetic which increases flexibility in 

combining operand types. The five standard operand types are as follows: 


COMPLEX 

DOUBLE 

REAL 

INTEGER 

LOGICAL 


two words per element 
two words per element 
one word per element 
one word per element 
one bit per element 


The programmer may define three non-standard types specifying multi-word 
elements or partial word elements, called bytes, whose length in bits is an integral 
divisor of 48. The mechanics of the TYPE declaration are covered in section 3.1. 


Mixed mode arithmetic is completely general; however, most applications will 
probably mix operand types real and integer, real and double, or real and complex. 
The following rules establish the relationship between the mode of an evaluated ex¬ 
pression and the types of the operands it contains. 


2 . 2.1 

TYPE AND 
MODE RULES 


AMI The order of dominance of the standard operand types within an expression from 
highest to lowest is: 

COMPLEX 

DOUBLE 

REAL 

INTEGER 

LOGICAL 


A M2 The mode of an evaluated arithmetic expression is referred to by the name 
of the dominant operand type. 

A M3 In mixed arithmetic expressions containing non-standard types the following 
restrictions hold: 

1. The non-standard types (types 5, 6, 7) may never be mixed with each other. 

2. Any one of the types 5, 6, 7 may be mixed with any or all of the standard 
types. When this is done, the non-standard type dominates the hierarchy 
established in rule AMI. 


AM4 In expressions of the form A**B the following rules apply: 

1. Neither A nor B may be type logical or byte (non-standard) type. 

2. B may be negative in which case the form is: A**(-B). 

3. For the standard types (except logical) the mode/type relationships are: 


Type A 

Type B 

I 

I 

R 

D 

C 

I 

R 

D 

C 

R 

R 

R 

D 

C 

D 

D 

D 

D 

C 

C 

C 

C 

C 

C 


mode of A**B 


For example, if A is real and B is complex, the mode of A**B is complex. 
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4. If A or B or both are a non-standard multi-word type, the programmer must 
provide subroutines for the evaluation of A**B. 

2 . 2.2 

MIXED MODE 

EVALUATION Examples 

1. Given A, B type real; I, J type integer. The mode of expression A*B-I+J 
will be real because the dominant operand is type real. It is evaluated: 

A*B—real 
Convert I to real 
R-l -I —*- R 2 real 
Convert J to real 

R 2 +J —►- R 3 real Evaluation completed 

2. The use of parentheses may change the evaluation. A,B,I,J are defined 
as above. A*B-(I-J) is evaluated: 

I-J-R x integer 

Convert R x to real-*~R 2 

A*B —R 3 real 

r 3 _r 2 Jkj real Evaluation completed 

3. Given C1,C2 type complex; A1,A2 type real. The mode of expression Al* 
(Cl/C2)+A2 will be complex because its dominant operand is type complex. 

It is evaluated: 

Cl/C2-*~Ri complex 

Convert Al to complex 

A1*R1-*~R2 complex 

Convert A2 to complex 

R 2 +A2 -*~R 3 complex Evaluation completed 

4. Consider the expression Cl/C2+(A1-A2) where the operands are defined as in 
3 above. It is evaluated: 

A1-A2 -*~Ri real 

Convert Rj to complex-»~R 2 
Cl/C2-*-R 3 complex 

Rj +R 2 -*-R 4 complex Evaluation completed 

Mixed mode arithmetic with standard types is illustrated by this example. 

5. Given: C complex 

D double 
R real 
I integer 
L logical 

and the expression C*D+R/l-L 
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The dominant operand type in this expression is type complex; therefore, the 
evaluated expression will be of mode complex. Evaluation: 


2.3 

MIXED MODE 
CONVERSIONS 


2.4 

ARITHMETIC 

REPLACEMENT 

STATEMENT 


Round D to a real and affix zero imaginary part 


C*D—*- Ri 
Convert R to 
R/I-^R 2 

Rs+Ri-Ra 

Convert L to 
R3 -L —R4 


complex 

complex; convert I to complex 

complex 

complex 

complex 

complex Evaluation completed 


If the same expression is rewritten with parentheses as C*D+(R/l-L) the 
evaluation proceeds: 


Convert I to 

real 

R/I--R 4 

real 

Convert L to 

real 

R 1 -L —*- R 2 

real 

Convert R 2 to 

complex—*-R 3 

Round D to 

real and affix zero imaginary part 

C*D —R 4 

complex 

R^ +R 3 Rg 

complex Evaluation completed 


Mixed mode arithmetic is accomplished through the special library subroutines. 

In the 1604 computer system, these routines include double precision and' complex 
arithmetic. In the 3600 system, the double precision arithmetic is built into the 
hardware; the complex arithmetic is performed by a library subroutine. 


The general form of the arithmetic replacement statement (or simply, arithmetic 
statement) is A = F, where F is an arithmetic expression and A is an identifier 
representing a variable. The operator = means that the value of the evaluated 
expression, F, is assigned to A, with conversion for mode if necessary. 

The identifier A is a variable; usually the type is a standard form: complex, double, 
real, integer, or logical. 

Non-standard types may also be specified; they may be used as left-hand variables 
also. 
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Complex and double precision variables are floating point quantities requiring 
two computer words. Real, integer and logical variables are represented by one 
word. The mode of an evaluated expression is determined by the type of dominant 
operand. However, this does not restrict the types that A may assume. An 
expression of complex mode may replace A even if A is of type real. The following 
chart shows the A,F relationship for all the standard modes. 


ARITHMETIC REPLACEMENT STATEMENT A = F 

A is an Identifier F is an Arithmetic Expression 
0(f) is the Evaluated Arithmetic Expression 


x. Mode of 

\ 0(f) 

TYPE x. 

of A N. 

Complex 

Double 

Real 

Integer 

Complex 

Store real & 
imaginary parts 
of 0(f) in 
real & imaginary 
parts of A. 

Round 0(f) to 
real. Store 
in real part of 

A. Store zero 
in imaginary 
part of A. 

Store 0(f) in real 
part of A. Store 
zero in imaginary 
part of A. 

Convert 0(f) 
to real & store 
in real part 
of A. Store 
zero in imaginary 
part of A. 

Double 

Discard imaginary 
part of 0(f) 

& replace it 
with ±0 according 
to real part of 0(f). 

Store 0(f) 

(most & least 
significant parts) 
in A (most & least 
significant parts). 

If 0(f) is ± affix 
±0 as least 
significant part. 
Store in A, most 
& least significant 
parts. 

Convert 0(f) 
to real. If 0(f) 
is ±, affix ±0 as 
least significant 
part. Store in A, 
most & least 
significant parts. 

Real 

Store real part 
of 0(f) in A. 
Imaginary part 
is lost. 

Round 0(f) to real 
& store in A. Least 
significant part of 

0(f) is lost. 

Store 0(f) in A. 

Convert 0(f) to 
real. Store in A. 

Integer 

Convert real part 
of 0(f) to INTEGER. 
Store in A. 

Imaginary part is 
lost. 

Round 0(f) to 
real, convert to 
INTEGER & store 
in A. The least 
significant part is 
lost. 

Convert 0(f) to 
INTEGER. Store 
in A. 

Store 0(f) in A. 

Logical 

If real part of 

0(f)^O, 1-^A. 

If real part of 

0(f)=O, 0-WY 

If 0(f)^O, 
store 1 in A. 

If 0(f)=O, store 

0 in A. 

Same as for 
double at left. 

Same as for double 
at left. 


When all of the operands in the expression F are of type logical, the expression 
is evaluated as if all the logical operands were integers. Let L x , L 2 , L 3 , L 4 
be logical variables, let R be a real variable and I an integer variable. 
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I — L 4 * L2+L3-L 4 


will be evaluated as if the L* were all integers (0 or 1) and the resulting value 
will be stored, as an integer, in I. 

R = L x *L 2 +L 3 -L 4 

is evaluated as stated above, but the result is converted to a real (a floating point 
quantity) before it is stored in R. 


Examples 



Given: Cj , 

Ai 

complex 

Di, 

a 2 

double 

Ri. 

A3 

real 

Ii , 

A4 

integer 

Li, 

As 

logical 

1. A, = C x 

*c 2 

-c 3 /c 4 


The mode of the expression is complex. Therefore the result of the expres¬ 
sion is a two-word, floating point quantity. A x is type complex and the result 
replaces A x . 

2. A 3 = Cx 

The mode of the expression is complex. The type of A 3 is real; therefore 
the real part of C 1 replaces \ . 

3. A 3 = 0^(0.,-l.) 

The mode of the expression is complex. The type of A 3 is real; the 
imaginary part of C 4 replaces A 3 . 

4. A4 =R 1 /R 2 *(R 3 -R 4 )+I 1 -(I 2 *R 5 ) 

The mode of the expression is real. The type of A 4 is integer; the result 
of the expression evaluation, a real, will be converted to an integer replacing 
A 4 . 

5. A 2 = d 4 **2 * (D 2 +(D 3 *D 4 ) )+(D 2 *R 1 *R 2 ) 

The mode of the expression is double. The type of A 2 is double; the result 
of the expression evaluation,, a double precision floating quantity replaces A 2 . 

6. A 5 = C^Rj-Ra+Ii 

The mode of the expression is complex. Since A 5 is type logical, an integer 
1 will replace A 5 if the real part of the evaluated expression is not zero. If 
the real part is zero, zero replaces A 5 . 


2.5 

LOGICAL AND 
RELATIONAL 

EXPRESSIONS A logical expression has the general form 

Ch op 0 2 op 0 3 ... op 

where Ch are arithmetic expressions, relations, or variables of type logical, 
and op is one of the logical operators .NOT. .AND. .OR. The value of a 
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logical expression is either true or false. A / relational expression has the form 
q iP q 2 where at least one of q 15 q 2 is an arithmetic expression; the other q may 
be either an arithmetic expression or a single logical variable, p is an operator 
belonging to the set 

.EQ. .NE. .GT. .GE. .LT. .LE. 

A relation is true if q x stands in the relation p to ^ . A relation is false if q 
does not stand in the relation p to . 

Within the compiler, relations are evaluated as illustrated in the following example. 
Consider the relation p = q. 

This is equivalent to the question, does p-q = 0? 

The compiler computes the difference and tests it for zero. If the difference is 
zero, the relation is true. If the difference is not zero, the relation is false. 

Relational expressions are converted internally to arithmetic expressions accord¬ 
ing to the rules of mixed mode arithmetic/*' These expressions are evaluated and 
compared with zero to determine the truth value of the corresponding relational 
expression. When expressions of mode complex are tested for zero, only the 
real part is used in the comparison. 


2.5.1 

RULES GOVERNING 

RELATIONS REL1 The only permissible forms of a relation are: 

q x P q 2 

q by itself, in which case a non-zero value is true and 

a zero value is false. 

REL2 q p q 2 p qg .. . is not permissible. 

REL3 The evaluation of a relation of the form q p q g is from left to right. The 
relations q x p q 2 , q x p (q^, (q i ) p q 2 , (q^ p (q g ) are equivalent. 


Examples 

A .GT. 16. R(I) .GE. R(I-l) 

R-Q(I)* Z .LE. 3.141592 K .LT. 16 

B-C .NE. D+E I .EQ. J(K) 

2.5.2 

LOGICAL 

EXPRESSION 

RULES LOG1 The hierarchy of logical operations is: 

First .NOT. 
then .AND. 
then .OR. 

LOG2 A logical variable or a relational expression is, in itself, a logical expression. 
If *C\, =zf 2 are logical expressions, then 

.NOT. JC x 

.AND. ^2 
-OR. c/C, 2 


17 




are logical expressions. If °C is a logical expression, (~/*), ( (^) ) are logical 
expressions. 

LOG3 If ^, 0^2 are logical expressions and op is .AND. or .OR. then, \ op op 
is never legitimate. 

L0G4 .NOT. may appear in combination with .AND. or .OR. only as follows: 

.AND..NOT. 

.OR.. NOT. 

.AND.(.NOT. • • •) 

.OR.(.NOT. • • •) 

.NOT. may appear with itself only in the form .NOT.(.NOT.(.NOT. • • ■ 

L0G5 If are logical expressions, the logical operators are defined as follows: 

.NOT. •=C 1 is false if is true 

.AND. is true if and only if are both true 

°C i .OR -c^ 2 ^ false if and only if a £\, B £ z are both false 


Examples 

Logical Expressions 


{The product A*B greater than 16.} .AND. (C equals 3.141519} 
A*B .GT. 16. .AND. C .EQ. 3.141519 
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{A(I) greater than 0} .OR. {B(J) less than 0 } 
A (I) .GT. 0 .OR. B(J) .LT. 0 



In the two examples below, all L t are of TYPE LOGICAL 
(L2 .OR. .NOT. L3) 



L2 .OR. .NOT, L3 .AND. (.NOT. L6 .OR. L5) 
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2.6 

LOGICAL 

REPLACEMENT 

STATEMENT 


2.7 

MASKING 

EXPRESSIONS 


2.7.1 

MASKING 

EXPRESSION 

RULES 


Incorrect Usages 

A.GT.(B.AND.C) 

10.LE.N.LE.100 
Q.NOT. .OR.R 
C.AND. .NOT. .NOT.B 

The last expression is permissible in the form C.AND. .NOT.(.NOT.B) 


The general form of the logical replacement statement is L=E, where L 
is a variable of type logical and E may be a logical expression or relation, or 
an arithmetic expression. 

When an arithmetic expression appears in a logical replacement statement, that 
expression is examined for being zero or non-zero. If the expression is non-zero, 
the left hand variable, L has the value TRUE. If the expression is equal to zero, 
the left hand variable L has the value FALSE. Thus the treatment of arithmetic 
expressions in logical replacement statements is consistent with that given to 
logical expressions in logical replacement statements. 


In FORTRAN-63, a masking expression is one in which 48-bit arithmetic is 
performed bit-by-bit on the operands within the expression. These operands 
must be of types real and integer only. Type integer includes octal and 
Hollerith constants. 

The masking operators are: .NOT. .AND. .OR. Although these operators are 
identical in appearance to the logical operators, their meanings are different. For 
masking operators the following definitions apply: 


.NOT. means complement the operand 

.AND. means form the bit-by-bit logical product of two operands 
.OR. means form the bit-by-bit logical sum of two operands 


The operations are described below. 


p 

V 

p .AND. v 

p .OR. v 

.NOT. p 

1 

1 

1 

1 

0 

1 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0 

1 


ME1 The hierarchy of operation is: first .NOT., then .AND., then .OR. 
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2.8 

MASKING 

REPLACEMENT 

STATEMENT 


ME2 Let B 1 ,B 2 be variables or constants whose types are real or integer. Then 
the following are masking expressions. 

.NOT. Bj 

B, .AND. B 
1 2 

B . OR. B 2 

ME3 If B is a masking expression, then (B), ((B) ),are masking expressions. 

ME4 .NOT. may appear with .AND. or .OR. only as follows: 

.AND..NOT. 

.OR..NOT. 

.AND. (.NOT. • • •) 

.OR. (.NOT. • • • ) 

ME5 Masking expressions of the following forms are evaluated from left to right. 

A .AND. B .AND. C . . . 

A .OR. B .OR. C . . . 


Examples 

Given: 


Ax 

7777000000000000 

octa 1 

A 2 

0000000077777777 

octa 1 

B 

0000000000001763 

octa 1 

C 

2004500000000000 

octa 1 


constant 

constant 

form of integer constant 
form of real 


.NOT. A 1 

A, .AND. C 

A y .AND. .NOT. C 
B .OR. .NOT. A 2 


is 0000777777777777 

is 2004000000000000 

is 5773000000000000 
is 7777777700001763 


The general form of the masking replacement statement is E = M. The masking 
statement is distinguished from the logical statement in the following ways. 


1. The type of E must be real or integer. 

2. All operands in the expression M must be type real or integer. M may 
contain parenthetical arithmetic subexpressions whose mode is real or 
integer. 
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Examples 

Given: All variables of type real or integer. 


2.9 

MULTIPLE 

REPLACEMENT 

STATEMENTS 


A(I) = B .OR. .NOT. C(I+2,J*K) 

B = D .AND. Q 

C(I,J) = .NOT. Z(K) .AND. (Q1 ..OR. .NOT. Q2) 
TEST = CELESTE .AND. 7HECLIPSE 
AB = D .OR. (S + T) 


The multiple replacement statement is a generalization of the replacement state¬ 
ments discussed earlier in this chapter, and its form is: 

'I'h = '/'n-i = . . . = ^ 2 = ^ 1 = expression 

The expression may be arithmetic, logical or masking. The 4 / l are variables 
subject to the following restrictions: 

Arithmetic or Logical Statement: 4 / 1 = EXP 

If EXP is logical or arithmetic, then 

If the variable ^ is type complex, double, real, or integer, 
then - EXP is an arithmetic statement. 

If the variable ^ is type logical, then \f/ l = EXP is a logical statement. 
Masking Statement: ^ = EXP 

If EXP is a masking expression, 4 / i must be a type real or integer variable 
only. 

The remaining n-1 ’/'i may be variables of any type, and the multiple replacement 
statement replaces each of the variables lA 2 , . . . , with the value of ^in a manner 

analogous to that employed in mixed mode arithmetic statements. 


Examples 


Given: A,B, C,D 

real 

The numbers in the examples 

E, F 

complex 

represent the evaluations of 

G,H 

double 

expressions. 

I,J 

integer 


K,L 

logical 
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I = A = 4.6 

A = I =4.6 

I = A = E = (10.2,3.0) 

F = A = I = E =(13.4,16.2) 

K = I = -14.6 

I = K = -14.6 


4.6—►- A 

4 —- I 

4 —- I 

4.0—►- A 

10.2—*- E real 
3.0 —*- E imaginary 

10.2- 

10 -- I 

13.4—E real 
16.2 —►- E imaginary 
13 —► I 

13.0—- A 
13.0—►- F real 
0.0—F imaginary 

-14 -*- I 

1 -K 

0 

1 -K 
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TYPE DECLARATIONS 
AND STORAGE ALLOCATION 


3 


3 


This chapter discusses how FORTRAN-63 allocates storage. The relation between 
word structure (TYPE) and array length (DIMENSION, COMMON) is explained. The 
methods for sharing storage (EQUIVALENCE) and the DATA statement is explained. 


3.1 

TYPE DECLARATIONS The TYPE declaration provides the compiler with information regarding the 

structure of the identifiers that name variables (1.6) and functions (5.1). The 
discussion that follows describes how type information is passed to the compiler 
from source language statements. 

There are five standard variable types (non-standard types are explained in 
Volume HI). Identifiers are declared of a given type by one of the following 
declarative statements: 


TYPE 

COMPLEX 

List 

TYPE 

DOUBLE 

List 

TYPE 

REAL 

List 

TYPE 

INTEGER 

List 

TYPE 

LOGICAL 

List 


A list, as used here, is a string of identifiers, in which each identifier is separated 
from the succeeding one by a comma. Subscripts are not permitted. An example 
of a list is: 

A,B1,CAT,D36F, EUPHORIA 
The characteristics of the standard variable types are: 


Type 

Element Definition 

Quantification 

Complex 

2 words/Element 

Floating point 

Double 

2 words/Element 

Floating point 

Real 

1 word /Element 

Floating point 

Integer 

1 word /Element 

Integer 

Logical 

1 bit /Element 

Logical 


3.1.1 

TYPE DECLARATION 

RULES TD1 The TYPE statement is non-executable and must precede the first executable state¬ 

ment in a given program. 

TD2 If a variable is declared differently in two or more TYPE statements, its TYPE 
will be determined from the last TYPE statement in which it appears. 
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3.2 

DIMENSION 


TD3 A variable not declared in a TYPE statement will be interpreted as TYPE REAL 

if the first letter of its identifier is A.H or O, . . .Z. It will be interpreted 

as TYPE INTEGER if the first letter of the identifier is I, J, K, L, M, N. 


Examples 

TYPE COMPLEX 
TYPE DOUBLE 
TYPE REAL 
TYPE INTEGER 
TYPE LOGICAL 


A147, RIGGISH, AT1LL2 
TEEPEE, B2BAZ 
EL, CAMINO, REAL, IDE63 
QUID, PRO, QUO 
GEORGE6 


A subscripted variable represents an element of an array of variables. Storage 
may be reserved for arrays by the non-executable statements DIMENSION or 
COMMON. 

The standard form of the DIMENSION statement is: 

DIMENSION Vi ,V 2 , . . . 

Vi have the form: Identifier (subscript string). The subscript string may have up 
to 3 unsigned constants separated by commas, as in SPACE(5,5,5). Under certain 
conditions within subprograms only, the subscripts may be integer variables. 
(Variable Dimensions 5.8) 

The number of computer words reserved for a given array is a function of the 
product of the subscripts in the subscript string, and the type of the variable. In 
the statements 

TYPE COMPLEX HERCULES 
DIMENSION HERCULES (10,20), 

the number of elements in the array HERCULES is 200. The TYPE statement, 
however, specifies two words per element; therefore, the number of computer words 
reserved is 400. The argument is the same for TYPE DOUBLE. For REAL and 
INTEGER the number of words in an array equals the number of elements in the 
array. 

For subscripted logical variables, up to 32 bits of a computer word are used; each 
bit represents an element of the logical variable array. The elements are stored 
left to right in a computer word starting with the most significant bit position. In 
the statements 

TYPE LOGICAL XERXES 
DIMENSION XERXES (5,5,5), 

there are 12 5 elements in the array XERXES and these elements will occupy four 
sequential words as shown below. 
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3.2.1 

VARIABLE DIMENSIONS 


3.3 

COMMON 


WORD 

WORD+1 

WORD+2 

WOKD+3 



Wlien an array identifier and its dimensions appear as formal parameters in a 
function or subroutine, the dimensions may be assigned through the actual param¬ 
eter list accompanying the function reference or subroutine call. The dimensions 
so assigned must not exceed the maximum array size specified by the DIMENSION 
statement in the calling program. See section 5.8, Variable Dimensions and 
Subprograms for details and examples. 


Just as an expression may contain sub-expressions, a program may contain, or 
call upon, subprograms (Chapter V). Such programs must be able to communicate, 
and they frequently require access to areas of information that they use in common. 
These areas are specified by the statement COMMON. The general form of this 
statement is: 

COMMON / (B) I x / List / (B) I 2 / List . . . 

B is a bank designator and has meaning only in the 3600 computer system. It is an 
unsigned integer constant between 0 and 7, and is ignored in a FORTRAN-63 program 
executed in the 1604 computer system. When B is omitted in 3600 programs, it is 
assumed to be the same as B = 0. 

I is a COMMON block identifier and it may be up to eight characters in length. It 
designates either labeled or numbered COMMON blocks and has the form: 

CiC 2 . . , C p 1< p< 8 

If Ci is alphabetic the identifier denotes a labeled COMMON block; the remaining 
characters may be alphabetic or numeric. If all the Ci are numeric, the identifier 
denotes a numbered COMMON block. If Ci is numeric, the remaining characters 
must be numeric. 


Examples 

Labeled COMMON Identifiers Numbered COMMON Identifiers 


AZ13 1 

MAXIMUS 146 

Z 3600 
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List has the form V x ,V 2 . . . , where Vj is of the form identifier (subscript string). 


COMMON A,B,C^ , 

COMMON/ /A,B,C,D T 

COMMON /BLOCKl/A,B/l234/C(10),D(10,10),E(10,10,10) 
COMMON/ (l)BLOCKA/D(l5) ,F(3,3),GOSH(2,3,4),Q1 


3.4 

COMMON BLOCKS The primary purpose of the COMMON block is to provide the programmer with a 

means of using, in subprograms, certain COMMON areas specified in the main 
program by referring only to the block desired. Both numbered and labeled blocks 
may be used for this purpose. Data stored in labeled COMMON blocks by the DATA 
statement are available to any subprogram using the appropriate labeled block. 


3.4.1 

COMMON RULES COM1 COMMON is non-executable and must precede the first executable statement in the 

program. 

COM2 If TYPE, DIMENSION or COMMON appear together, the order is immaterial. 

COM3 The identifiers of labeled COMMON blocks are used only for block identification 
within the compiler; they may be used elsewhere in the program as other kinds of 
identifiers. 

COM4 For any given dimensioned variable, the dimensions may be declared either in a 
COMMON statement or in a DIMENSION statement. If declared in both, those of 
the DIMENSION statement override those declared hi the COMMON statement. 

COM5 At the beginning of program execution, the contents of the COMMON area are 
undefined unless specified by a DATA statement. 

COM6 An identifier in one COMMON block may not appear in another COMMON block. If 
it does, the identifier is' doubly defined. 


3.4.2 

COMMON BLOCK LENGTH The length of a COMMON block, in computer words, is determined from the type 

of the list identifier and the dimension (if any) associated with that identifier. 

Given COMMON/A/Q(4),R(4),S(2) 

TYPE COMPLEX S 

the length of the COMMON block A is 12 computer words. The origin of the 
COMMON block is Q(l). 


t 


These forms are sometimes calle 


d blank COMMON. 
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block A 


origin 


Q(l) 

Q(2) 

Q(3) 

Q(4) 

R(l) 

R(2) 

R(3) 

R(4) 

S(l) 

S(l) 

S(2) 

S(2) 



real part 
imaginary part 
real part 
imaginary part 


Examples 

MAIN PROG 


TYPE COMPLEX C 
COMMON/TEST/C(20)/36/A,B,Z 


The length of TEST is 40 computer words. 


The subprogram may re-arrange the allocation of words as in: 


SUBPROG1 


COMMON/TEST/A(10),G(10),K(10) 
TYPE COMPLEX A 


The length of TEST is 40 words. The first 10 elements (20 words) of the block, 
represented by A, are complex elements. Array G is the next 10 words, and array 
K is the last 10 words. Within the subprogram, elements of G will be treated as 
floating point quantities; elements of K will be treated as integer quantities. 


The length of the COMMON block must not be changed by the subprograms using the 
block. The identifiers used within the block may differ as shown above. 


The following arrangements are equivalent: 


( TYPE DOUBLE A 

< DIMENSION A(10) 
( COMMON A 

( DIMENSION A(10) 

< TYPE DOUBLE A 
t COMMON A 

( COMMON A 

< DIMENSION A(10) 
l TYPE DOUBLE A 


( TYPE DOUBLE A 
< COMMON A 
l DIMENSION A(10) 

( TYPE DOUBLE A 
l COMMON A(10) 


The label of a COMMON block is used only for block identification. The following 
is not erroneous: 


COMMON /A/A(10)/B/B(5j5) /C/C (5,5,5) 







3.5 


EQUIVALENCE 


3.5.1 

EQUIVALENCE RULES 


The EQUIVALENCE statement permits variables to share locations in storage. 
The general form of this statement is: 

EQUIVALENCE (A,B . . .), (A1,B1, 

where the A, B, ... are simple or singly subscripted variable identifiers. A 
multiply subscripted variable can be represented by a singly subscripted variable. 
The correspondence is: 

A (i,j,k) = A (the value of (i+(j-l)*I + (k-l) *I*J) ) 

where i,j,k are integer constants and I and J are the integer constants appearing 
in DIMENSION A (I,J,K). For example, given DIMENSION A(2,3,4), the element 
A(l,l,2) is represented by A(7). 


EQU1 EQUIVALENCE is non-executable and must precede the first executable statement 
in the program or subprogram in which it appears. 

EQU2 If TYPE, DIMENSION, COMMON, or EQUIVALENCE appear together, the order 
is immaterial. 


EQU3 The following may be made equivalent 

COMPLEX / COMPLEX 

COMPLEX / DOUBLE 

COMPLEX or DOUBLE / REAL 

COMPLEX or DOUBLE / INTEGER 

REAL / REAL 

REAL / INTEGER 

DOUBLE / DOUBLE 

INTEGER / INTEGER 

LOGICAL / LOGICAL 

TYPE 5 / TYPE 5 (non-standard) 

TYPE 6 / TYPE 6 (non-standard) 

TYPE 7 / TYPE 7 (non-standard) 


with or without 
subscript 


unsubscripted only 


♦ 


Any variable of TYPE LOGICAL, 5, 6, or 7 may be made equivalent to one of the 
standard types, but they must not be subscripted. 

EQU4 The EQUIVALENCE statement does not rearrange COMMON, but arrays may be 
defined as equivalent so that the length of the COMMON block is changed. The 
origin of the COMMON block must not be changed by the EQUIVALENCE statement. 

The following simple cases illustrate changes in block lengths caused by the 
EQUIVALENCE statement. 


Given: Arrays A and B 

S a = subscript of A 
St, = subscript of B 
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CASE I A, B both in COMMON 



a) 


b) 


If A appears before B in the COMMON statement: 

Sa — Sb is a permissible subscript arrangement 
Sa < Sb is not 

If B appears before A in the COMMON statement 

Sa 2 = sb is a permissible subscript arrangement 
Sa > Sb is not 

Block 1 


A 

(1) 



A 

(2) 

B 

(1) 

A 

CO 

B 

(2) 

A 

(4) 

B 

(3) 

A 

(5) 

B 

(4) 



B 

(5) 



B 

(6) 



B 

(7) 


COMMON/l/A(5), B (7) 
EQUIVALENCE (A(4), B(3) ) 



Statement EQUIVALENCE (A(3), B(4) ) changes the origin of block 1. This is not 
permitted. 

B(l) —origin changed 
origin—^A(l) B(2) 

A(2) B(3) 

A(3) B(4) 

A(4) B(5) 


CASE II A in COMMON, B not in COMMON (corresponds to CASE la) 


Sb £ Sa is a permissible subscript arrangement 


Sb > Sa is not 


Block 1 

A(l) 


COMMON /l/A(4) 

A(2) 

B(l) 

DIMENSION B(5) 

A (3) 

B(2) 

EQUIVALENCE (A(3), B(2) 

A(4) 

B(3) 



B(4) 

B(5) 



CASE III B in COMMON, A not in COMMON (corresponds to CASE lb) 


Sa ^ Sb is a permissible subscript 
Sa > Sb is not 



Block 1 

origin—►- B(l) 


COMMON/1/ B (4) 

B(2) 

A(l) 

DIMENSION A (5) 

B(3) 

A(2) 

EQUIVALENCE (B(2), A(l) 

B(4) 

A(3) 



A(4) 

A(5) 



CASE IV A, B, not in COMMON 

No subscript arrangement restrictions. 
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3.5.2 

GENERAL 


3:6 

DATA 


RULE Regarding EQUIVALENCE and COMMON - Consider the statement 

EQUIVALENCE (A(6),B(4),C(3) ,D(8) ) 

The base of the equivalence is the identifier with the largest subscript. The base 
is D(8); A(6),B(4), and C(3) will be made equivalent to it. 

If any, or all, of A, B, C, D occur in a COMMON statement, the order, from left to 
right, is by descending subscripts in the EQUIVALENCE statement. Since the 
subscript of D is greater than the subscript of A, et cetera, the following COMMON 
statement is permissible: 


COMMON/l/D(10),A(8),B(5),C(10) 


The combined statements 

EQUIVALENCE (A(6),B(4),C(3),D(8) ) 
COMMON/l/D(10),A(8),B(5),C(10) 


yield the storage arrangement: 

Block 1 


origin 

D(l) 





D(2) 





D(3) 

A(l) 




D(4) 

A (2) 




D(5) 

A(3) 

B(l) 



D(6) 

A(4) 

B(2) 

C(l) 


D(7) 

A(5) 

B(3) 

C(2) 

equivalence base 

D(8) 

A(6) 

B(4) 

C(3) 


D(9) 

A(7) 

B(5) 

C(4) 


D(10) 

A(8) 


C(5) 





C(6) 





C(7) 





C(8) 





C(9) 





C(10) 


Within the EQUIVALENCE statement, the order is immaterial. EQUIVALENCE 
(A(6),B(4),C(3),D(8) ) 

is the same as EQUIVALENCE(A(6),D(8),C(3) ,B(4) ). 


The programmer may assign constant values to variables in the program by using 
the DATA statement either by itself or with a DIMENSION statement. It may be 
used to store constant values in variables contained in a labeled COMMON block. 

The form of the DATA statement is: 

DATA(Ii =List),(I 2 =List), ... 

List contains constants only and has the form 

a x ,a 2 , . . . ,K(bj, ,b 2 , . . ,c 2 , . . . 

where K is an integer constant repetition factor that causes the parenthetical list 
following it to be repeated K times. I is an identifier representing a simple 
variable, a variable with integer constant subscripts, an array, or an array with 
integer variable subscripts. 
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3.6.1 

DATA RULES 


DAT1 DATA is non-executable and must precede the first executable statement in any 
program or subprogram in which it appears. 

DAT2 When DATA appears with TYPE, DIMENSION, COMMON or EQUIVALENCE state¬ 
ments, the order is immaterial. 

DAT3 DO loop-implying notation is permissible with the restriction that m 3 cannot appear. 
Short notation may be used for storing constant values in arrays. 

DAT4 No array name declared in blank or numbered COMMON or in a variable DIMENSION 
can belong to a DATA statement. 

DAT5 When a signed constant appears in a DATA list, that sign is unary; in DATA 

(A = -2.), the negative value of the floating point number 2 replaces A. Negative 
octal constants are prefixed with minus signs. The operator.NOT.may not 
be used. 

DAT6 With identifiers of types real or integer, the corresponding constant in the list must 
be the same type; in DATA (A = 2). The type of A is not checked, and an integer 2 
will replace A. 

DAT7 There must be a one-one correspondence between the identifier and the list. This 
is particularly important in arrays. 

Consider COMMON /BLK/ A(3), B 

DATA (A = l.,2.,3.,4.) 

The constants l.,2.,3., are stored in array locations A, A+l, A+2; the constant 4, is 
stored in location B. If this occurs unintentionally, erroneous results may occur 
when B is referred to elsewhere in the program. 

Consider COMMON / TUP/ C(3) 

DATA (C = 1..2.) 

The constants 1., 2. are stored in array locations C and C+l, the contents of C(3), 
that is, location C+2 are not defined. 

DATS Use of DATA with a TYPE LOGICAL variable constitutes a special case (the last 
example below). 


Examples 

DATA (LEDA=15), (CASTOR=16.0), (POLLUX=84.0) 
LEDA 15 


CASTOR 16.0 

POLLUX 84.0 

DATA (A(l,3) = 16.239) 

ARRAY A 

A(l,3) 16.239 
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DIMENSION B(10) 

DATA (B = 77B, -77B, 4(776B, -774B) ) 

ARRAY B 77B 

-77B 
77 6B 
-774B 
77 6B 
-774B 
77 6B 
-774B 
776B 
-774B 

COMMON /HERA/ C(4) 

DATA (C = 3.6, 3(10.5) ) 

ARRAY C 3.6 

10.5 

10.5 

10.5 

TYPE COMPLEX PROTEUS 

DIMENSION PROTEUS (4) 

DATA (PROTEUS = 4( (1.0, 2.0) ) ) 

ARRAY PROTEUS 1.0 

2.0 

1.0 

2.0 

1.0 

2.0 

1.0 

2.0 

DIMENSION MESSAGE (3) 

DATA (MESSAGE = 3HWHO, 2HIS, 6HSYLVIA) 

ARRAY MESSAGE WHO 

IS 

SYLVIA 

This example illustrates how elements of a logical array are stored by the DATA 

statement. 

Given: TYPE LOGICAL L 

COMMON / NETWORK / L (4,8) 

Store the following matrix of logical elements: 

11101011 
10110100 
L “ 11110011 
10100010 
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Arrays are stored columnwise. 

Elements of logical arrays are stored 32 bits to the word, left to right, left jus¬ 
tified with zero fill. 

The matrix fits into one computer word as follows: 

111 110 101 111 011 010 000 100 101 110 100 0. . . 0 
and its octal equivalent is 
7657320456400000 

Therefore, the appropriate DATA statement is: 

DATA (L = 7657320456400000B) 
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4.1 

STATEMENT 

IDENTIFIERS 


4.2 

GO TO STATEMENTS 


4.2.1 

UNCONDITIONAL GO TO 


4.2.2 

ASSIGNED GO TO 


CONTROL STATEMENTS 



Program execution normally proceeds from one statement to the statement 
immediately following it in the program. Control statements can be used to alter 
this sequence or cause a number of iterations of a program section. 

Control may be transferred to an executable statement only; a transfer to a 
non-executable statement will result in a program error. During compilation, 
however, no error will be indicated. 

Iteration control provided by the DO statement causes a predetermined sequence 
of instructions to be repeated any number of times with the stepping of a simple 
integer variable after each iteration. 


Statements are identified by numbers which can be referred to from other 
sections of the program. A statement number used as a label or tag, appears 
in columns 1 through 5 on the same line as the statement on the coding form. 
The statement number N may lie in the range 1 ^N^99999. An N of fewer than 
5 digits may occupy any of the first five columns; blanks are squeezed out and 
leading zeros are ignored, 1=01=001=0001, (Appendix A). 


Unconditional transfer of control is provided by GO TO statements. 


GO TO n 

This statement causes an unconditional transfer to the statement labeled n; 
n is a statement number. 


GO TO m, (n x ,n 2 , . . n m ) 

This statement acts as a many-branch GO TO. 

m is an integer variable assigned an integer value iq in a preceding ASSIGN 
Statement. The n 4 are statement numbers. The parenthetical list need not 
be present. 
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4.2.3 

ASSIGN STATEMENT 


4.2.4 

COMPUTED GO TO 


4.3 

IF STATEMENTS 


4.3.1 

THREE BRANCH IF 
(ARITHMETIC) 


ASSIGN A. TO in 

This statement is used with the Assigned GO TO statement, 
a statement number, m is a simple integer variable. 

ASSIGN 10 TO LSWTCH 


GO TO LSWTCH,(5,10,15,20) 
Control would transfer to statement 10. 


GO TO (n 1; n 2 , . . ., n m ), i 

This statement acts as a many-branch GO TO where i is preset or computed 
prior to its use in the GO TO. 

The iij are statement numbers and i is a simple integer variable. If i— 1, a 
transfer to n 1 occurs; if i— m, a transfer to occurs. 

ISWITCH = 1 

GO TO (10,20,30),ISWITCH 


10 JSWITCH = ISWITCH +1 
GO TO (11,21,31),JSWITCH 

Control would transfer to statement 21. 


0]s 

Conditional transfer of control is provided by the two- and three-branch IF 
statements, the status of sense lights or switches, or the status of an arithmetic 
overflow indicator. 


IF (A) n^ngpig 

A is an arithmetic expression and the n^ are statement numbers. 

This statement tests the evaluated quantity A and jumps according to the 
following criteria: 

A < 0 jump to statement n x 

A = 0 jump to statement n 2 

A > 0 jump to statement n 3 

In the test for zero, +0 = -0. When the mode of the evaluated expression is 

complex, only the real part is tested for zero. 

IF(A*B-C*SINF(X) )10,10,20 

IF(I)5,6,7 

IF(A/B**2)3,6,6 


V 
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4.3.2 

TWO BRANCH IF 
(LOGICAL) 


4.3.3 

SENSE LIGHT 


4.3.4 

SENSE SWITCH 


4.4 

FAULT CONDITION 
STATEMENTS 


IF(L) 

'L is a logical or an arithmetic expression. The iij are statement numbers. 
The evaluated expression is tested for true (non-zero) or false (zero). If L is 
true jump to statement n 1 . If L is false jump to statement n 2 . 

IF (A jGT. 16. .OR. I .EQ. 0)5,10 

IF(L)1,2 (L is TYPE LOGICAL) 

IF(A*B-C)1,2 (A*B-C is arithmetic) 

IF(A*B/C .LE. 14.32)4,6 

In the statement IF (A) 2,3,4, 

A is tested as shown in 4.3.1. In the statement IF (A) 4,3 

if A is not zero, jump to statement 4; if A is zero, jump to statement 3. 


IF (SENSE LIGHT i)n y ,n 2 

The statement tests sense light i. If it is on, it is turned off and a jump occurs 
to statement n r If it is off, a jump occurs to statement n 2 . 

i is a sense light and the n t are statement ntimbers. i may be a simple integer 
variable or constant. 

IF(SENSE LIGHT 4)10,20 


IF (SENSE SWITCH i)n x ,n 2 

If sense switch i is set (ON) a jump occurs to statement . If it is not set (OFF) 
a jump occurs to statement n 2 ; i may be a simple integer variable or constant. 

In the 3600 1 < i < 6 (physical console switches) 

In the 1604 1 < i < 6 (CO OP Monitor function. Appendix E) 

N = 5 

IF(SENSE SWITCH N)5,10 


At execute time the computer is set to interrupt on divide, overflow or exponent 
fault. 

IF DIVIDE CHECK nj,n 2 
IF DIVIDE FAULT n^ 

A divide fault occurs following division by zero. The statement checks for this 
fault; if it has occurred, the indicator is turned off and a jump to statement n 1 takes 
place. If no fault exists, a jump to statement n 2 takes place. 

IF EXPONENT FAULT n l5 n 2 

An exponent fault occurs when the result of a real or double or complex arithmetic 
operation exceeds the upper limits specified for these types. Results that are less 
than the lower limits are set to zero without indication. This statement is therefore 
a test for floating-point overflow only. If the fault has occurred, the indicator is 
turned off, and a jump to statement n 1 takes place. If no fault exists a jump to 
statement n 2 takes place. 
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IF OVERFLOW FAULT n^ 


4.5 

DO STATEMENT 


4.5.1 

DO INDEX VARIABLE: i 


An overflow fault occurs when the magnitude of the result of an integer sum or 
difference exceeds 2 47 -1. This fault does not occur in division and it is not 
indicated in multiplication. If the fault occurs, the indicator is turned off and a jump 
to statement n x takes place. If no fault exists, a jump to statement n 2 takes place. 


DO n i = m 1 ,m 3 

The DO Statement provides FORTRAN-63 with a recursive property. 

n is a statement number; i is the index variable. It is a simple integer variable, 
mi are the indexing parameters; they may be unsigned integer constants or simple 
integer variables. If m 3 does not appear, it is construed to be 1. 

The DO Statement, the statement labeled n, and any intermediate statements 
constitute a DO loop. Statement n may not be an IF or GO TO statement or 
another DO statement. The statement immediately following the DO statement 
must be executable (Appendix C). 


The initial value of i is m x . This value is compared with m 2 before executing the 
DO loop and, if it does not exceed m 2 , the loop is executed. After this step, i is 
increased by m 3 and control passes to the top of the loop where i is again compared 
with m 2 ; this process continues until i exceeds m 2 as shown below. Control then 
passes to the statement immediately following n, and the DO loop is said to be 
satisfied. Should m x exceed m 2 on the initial entry to the loop, the loop is not 
executed and control passes to the next statement. 



When the DO loop is satisfied, the index variable i is no longer well defined. If a 
transfer out of the DO loop occurs before the DO is satisfied, the value of i is 
preserved and may be used in subsequent statements. 
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DO NESTS 


When a DO loop contains another DO loop the grouping is called a DO nest. If 
Di ,D 2 , . . ,D m represent DO statements, where the subscripts indicate that D x appears 
before D 2 appears before D 3 , et cetera, and n x ,ng , , . . , 1 ^ represent the corresponding 
limits of the Di, then n^ must appear before n m _ ] ... n 2 must appear before n t . 


DO loops may be nested in common with other DO loops: 



Examples 

DO 1 1= 1,10,2 


DO 100 L=2,LIMIT 


DO 5 1=1,5 
DO 5 J=I,10 
DO 5 K=J,15 


DO 2 J=l,5 


DO 10 1=1,10 
DO 10 J=l,10 


5 CONTINUE 


DO 3 K=2,8 


3 CONTINUE 


2 CONTINUE 


10 CONTINUE 


DO 20 K=K1,K2 


20 CONTINUE 


DO 4 L=l,3 


100 CONTINUE 


4 CONTINUE 


1 CONTINUE 
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DO LOOP TRANSFER 


4.5.4 

DO PROPERTIES 


4.6 

CONTINUE 


In a DO nest, a transfer may be made from one DO loop into a DO loop that contains 
it; and a transfer out of a DO nest is permissible. The special case is transferring 
out of a nested DO loop and then transferring back to the nest. 

In a DO nest: 

If the range of i includes the range of j and a transfer out of the range of j occurs, 
then a transfer into the range of i or j is permissible. 

In the following diagram, EXTR represents a portion of the program outside of the 
DO nest. EXTR must not change the indexing variable or the indexing parameters. 



1) The indexing parameters m x , m 2 , m 3 are either integer constants or simple 
integer variables. 

2) The values of the indexing parameters are assumed to remain constant until the 
DO is satisfied. 

3) The indexing parameters should assume positive values only. 

4) If m x > lx^ initially, the loop is not executed. 

5) The identity and value of the indexing variable is local to the statements in the 
range of the DO statement when 

(a) it is not used as an operand 

(b) No transfers out of the range of the DO exist 
Otherwise, the identity and value of i is global. 

6) DO-loops may be nested 50 deep. 


CONTINUE 

The CONTINUE statement is most frequently used as the last statement of a DO loop 
to provide a transfer address for IF and GO TO instructions that are intended to 
begin another repetition of the loop. If CONTINUE is used elsewhere in the source 
program, it acts as a do-nothing instruction and control passes to the next sequential 
program statement. 



4.7 

PAUSE 


PAUSE 



4.8 

STOP 


4.9 

END 


PAUSE n 

n is an octal number such that 1 < n < 2 47 -1. PAUSE n halts the computer with n 
displayed in the accumulator register on the console. When the START key on the 
console is pressed, program execution proceeds starting with the statement 
immediately following PAUSE. 


STOP 
STOP n 

n is an octal number such that 1 <n < 2 47 -1. STOP n halts the computer with n 
in the accumulator register displayed on the console. When the START key on the 
console is pressed, an exit will be made to the COOP MONITOR (1604) or SCOPE 
(3600). STOP (n omitted) causes immediate exit to MONITOR or SCOPE. 


END 

END marks the physical end of a program or subprogram; if executed, it acts as a 
RETURN. 




43 




FUNCTIONS AND SUBROUTINES 


5 


5.1 

STATEMENT 

FUNCTIONS 


5.1.1 

STATEMENT 
FUNCTION RULES 


FORTRAN-63 functions and subroutines, range from single source language state¬ 
ments to independently compilable subprograms. 

A function name is constructed in the same way as a variable identifier and has a 
type determined by the conventions established for variables. A function together 
with its arguments may be used at any place in an expression that a variable 
identifier may be used. 

A reference to a function is a call upon a computational procedure for the return of 
a single value, identified by and associated jvith the function identifier. This 
procedure may be defined in a single statement within the program (statement 
function); it may be defined within the compiler (library function) or it may be 
defined in a multi-statement subprogram either compiled with a main program or 
compiled independently (function subroutine). 

A reference to a subroutine is also a call upon a computation procedure. This 
procedure may return one or more values or it may return none. No value is 
associated with the name of the subroutine, and the subroutine must be called by a 
CALL statement. 

Any function reference must supply the function with a set of arguments or parameters. 
This set must contain at least one argument and may contain up to 63 arguments. 

The forms of the arguments differ somewhat in each of the three kinds of functions. 

The form of the function reference is: 

F (p 15 p 2 , . . .p n ) 1 -n ^63 

where F is the function name and the pi are function arguments or actual parameters. 
The corresponding arguments appearing with the function name in a function 
definition are called formal parameters. 


Statement functions are defined by a single arithmetic or logical statement in the 
source program and apply only to the particular program or subprogram in which 
the definition appears. They have the form 

F (p 1 ,p 2 ,- . .p n ) = E 

where F is the function name and E is an expression. 


SF1 The type of the function is determined from the naming conventions specified for 
variables in Chapter 3, Type Declarations. 

SF2 The function name must not appear in a DIMENSION, EQUIVALENCE or COMMON 
statement. 
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SF3 The formal parameters will usually appear in the expression E. When the state¬ 
ment function is executed, the formal parameters are replaced by the corresponding 
act ual parameters of the function reference. Each of the formal parameters may be 
TYPE REAL or INTEGER only, but they may not be declared in a TYPE statement. 
(3.1.1, rule TD3) Each of the actual parameters may be any arithmetic expression, 
but there must be agreement in order, number and type between the actual and 
formal parameters. 

SF4 E may be arithmetic or logical. 

SF5 E cannot contain subscripted variables. 

SF6 The expression E may refer to library functions, previously defined statement 
functions and function subprograms. 

SF7 All statement functions must precede the first executable statement of the program 
or subprogram, but they must follow all declarative statements (DIMENSION, TYPE, 
et cetera). 


Examples 

TYPE COMPLEX Z 

Z (X,Y)=(1. ,0.) *EXPF(X) *COSF(Y)+(0. ,1.) *EXPF(X) *SINF (Y) 

This arithmetic statement function computes the complex exponential Z(x,y) = 
e x+iy . 


5.2 

LIBRARY FUNCTIONS FORTRAN-63 contains the standard library functions available in earlier versions 

of FORTRAN. A list of these functions is in Appendix D. The identifying names 
have not been changed. When one appears in the source program, a special part of 
the compiler identifies it as a library function and takes appropriate action as 
explained below. 

In Chapter 3, specific rules are given for declaring the types of identifiers. In the 
absence of a TYPE declaration, a variable type is determined by its first identifier 
letter. As stated in rule SF1, this convention applies to function identifiers. In the 
standard library function for obtaining the natural logarithm of a number (LOGF) 
the first identifier letter, L, would cause that function to return an integer result. 
In this case the result is contrary to established FORTRAN usage. To avoid 
inconsistency, the compiler recognizes the standard library functions and permits 
the programmer to use such functions in the usual manner. 

5.3 

FUNCTION 

SUBPROGRAMS Function subprograms are FORTRAN source language programs that cannot be 

defined by one statement and that are not used frequently enough to be included in 
the library. 

Function subprograms may be compiled independently, and the first statement of 
such a subprogram must have the form: 

FUNCTION F (p 15 p 2 ,. • -P n ) 1 —n — 63 
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where F is the function name, and the pj are formal parameters. These parameters 
may be array names, non-subscripted variables, or names of other function or 
subroutine subprograms. 


5.3.1 

FUNCTION SUBPROGRAM 

RULES FS1 The type of the function is determined from the naming conventions specified for 

variables in Chapter 3, Type Declarations. 

FS2 The name of a function must not appear in a DIMENSION statement. The name 
must appear, however, at least once as any of the following: 

The left-hand identifier of a replacement statement 

An element of an input list 

An actual parameter of a subprogram call 

FS3 No element of a formal parameter list may appear in a COMMON or EQUIVALENCE 
statement within the function subprogram. „ 

FS3A When a formal parameter represents an array, it should be declared in a 

DIMENSION statement within the function subprogram. If it is not declared, only 
the first element of the array will be available to the function subprogram. 

FS4 In referring to a function subprogram the following forms of the actual parameters 
are permissible: 

a. arithmetic expression 

b. constant or variable, simple or subscripted 

c. array name 

d. function reference 

e. subroutine 

In form d, if only the name of the function appears, it must also appear in an 
EXTERNAL statement in the calling program. See example 2 following 5.5. 

In form e, the subroutine may appear as a subroutine name alone or as a subroutine 
name with a parameter bit: 

These cases are illustrated in examples 3 and 4 following 5.5. 

FS4A Logical expressions may not be actual parameters. 

FS5 Actual and formal parameters must agree in order, number and type. 


FUNCTION TYPE 

AND MODE The compiler distinguishes between array names and functions as follows: Given 

an identifier followed by a left parenthesis, Z( , if the identifier Z occurs in a 
DIMENSION statement, it represents an array. If not, Z( represents a function. 
To determine the mode of the function (statement, subprogram, or library), the 
compiler goes through the following process of elimination: 

1. If the identifier is not in the compiler's table of library functions, or is 
not declared in a TYPE statement, the mode of the evaluated function is 
according to the identifier first-letter criterion (3.1.1). 
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5.4 

RETURN AND 
END STATEMENTS 


5.5 

EXTERNAL 

STATEMENT 


2. If the identifier is in the compiler's table of library functions, but not in 
a TYPE statement, the mode of the evaluated function is given by the 
library function (Appendix D). 

3. If the name is declared in a TYPE statement, the mode of the evaluated 
function is defined by the TYPE statement (3.1.1). 


A subprogram normally contains one or several RETURN statements that indicate 
the end of logic flow within the subprogram, and return control to the calling pro¬ 
gram. In function references, control returns to the statement in which the function 
is imbedded. In subroutine subprograms, control returns to the next executable 
statement immediately following the CALL statement in the calling program. The 
form of this statement is RETURN. 

The END statement marks the physical end of a program, subroutine subprogram 
or function subprogram. If the RETURN statement is omitted, END acts as a 
return to the calling program. 


When the actual parameter list of a given function reference contains a function or 
subroutine name, that name must be declared in an EXTERNAL statement. Its 
form is: 

EXTERNAL identifier, identifier.,, . . . 

where identifier is the name of a function or subroutine. The EXTERNAL statement 
must precede the first executable statement of any program in which it appears. 
When it is used, EXTERNAL always appears in the calling program. 


Examples 

1. Function Subprogram 

FUNCTION GREATER (A,B) 

IF (A .GT. B) 1,2 

1 GREATER = A-B 
RETURN 

2 GREATER = A+B 
END 

Calling Program Reference 

Z(I,J) = F1+F2-GREATER (C-D,3.*I/J) 

2. Function Subprogram 

FUNCTION PHI (ALFA, PHI2) 

PHI = PHI2(ALFA) 

END 
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Calling Program Reference 


EXTERNAL SINF 


C=D-PHI (Q(K) ,SINF) 

Frona its call in the main program, the formal parameter ALFA is 
replaced by Q(K), and the formal parameter PHI2 is replaced by SINF. 
PHI will be replaced by the sine of Q(K). 

3. Function Subprogram 

FUNCTION PSYCHE (A,B,X) 

CALL X 

PSYCHE = A/B*2*(A-B) 

END 

Function Subprogram Reference 
EXTERNAL EROS 


R = S - PSYCHE (TLIM,ULIM,EROS) 

In the function subprogram, TLIM, ULIM replaces A,B. The CALL X is 
a call to a subroutine named EROS. EROS appears in an EXTERNAL 
statement so that the compiler recognizes it as a subroutine name rather 
than a variable identifier. 

4. Function Subprogram 

FUNCTION AL(W,X,Y,Z) 

CALL W(X,Y,Z) 

AL = Z**4 
RETURN 

Function Subprogram Reference 
EXTERNAL SUM 


G = AL(SUM,E,V,H) 

In the function subprogram the name of the subroutine (SUM) and its 
parameters (E,V,H) replace W and X,Y,Z. SUM appears in the EXTERNAL 
statement so that the compiler will treat it as a subroutine name rather 
than a variable identifier. 
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5.6 


SUBROUTINE 

SUBPROGRAMS Subroutine subprograms may be compiled independently; the first statement of 

such a program must have the form: 

SUBROUTINE S 
or 

SUBROUTINE S ( Pl ,p 2 , . . ,p n ) 1 ^n ^63 

where S is the subroutine name, and the Pi are the formal parameters which may 
be array names, non-subscripted variables, or names of other function or sub¬ 
routine subprograms. 


5.6.1 

SUBROUTINE RULES SSI The name of the subroutine may not appear in any declarative statement (TYPE, 

DIMENSION) in the subroutine. 

552 The name of the subroutine must never appear within the subroutine as an 
identifier in a replacement statement, in an input-output list, or as an argument 
of another CALL. 

553 No element of a formal parameter list may appear in a COMMON or EQUIVALENCE 
statement within the subroutine subprogram. 

SS3A When a formal parameter represents an array, it should be declared in a 

DIMENSION statement within the subroutine. If it is not declared, only the first 
element of the array will be available to the subroutine. 

The executable statement in the calling program for referring to a subroutine 
subprogram is of the form: 

CALL S 
or 

CALL S (p,^, - * . Pn ) 1 — n — 63 

where S is the subprogram name, and the Pi are the actual parameters. 


5.6.2 

SUBROUTINE REFERENCE 

RULES SS4 The subroutine returns values through parameters or COMMON variables. No 

value is associated with its name. 

555 The subroutine name may not appear in any declarative statement (TYPE, 
DIMENSION et cetera). 

556 In the subroutine call, the following forms of actual parameters are permissible: 

a. arithmetic expression 

b. constant or variable, simple or subscripted 

c. array name 

d. function reference 

e. subroutine 

In form d, if only the name of the function appears, it must also appear in an 
EXTERNAL statement in the calling program (example 2, 5.5). In form e. the 
subroutine may appear as a subroutine name alone or as a subroutine name with 
a parameter list. 
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SS6A Logical expressions may not be actual parameters. 

SS7 Actual and formal parameters must agree in order, number and type. 

Examples 

1. Subroutine Subprogram 

SUBROUTINE BLVDLDR (A,B,W) 

W - 2. *B/A 
END 

Calling Program References 

CALL BLVDLDR (X(I),Y(I),W) 


CALL BLVDLDR (X(I)+H/2.,Y(I)+C(1)/2.,W) 


CALL BLVDLDR (X(I)+H,Y(I)+C(3) ,Z) 

2. Subroutine Subprogram (Matrix Multiply) 

SUBROUTINE MATMULT 
COMMON/BLK1/X(20,20),Y(20,20),Z(20,20) 
DO 10 1=1,20 

DO 10 J=l,20 

Z(I,J) = 0. 

DO 10 K=l,20 

10 Z(I,J) = Z(I,J) +X(I,K) *Y (K,J) 

RETURN 

END 

Calling Program Reference 

COMMON/BLKl/A (2 0,20), B (20,20), C (20,20) 


CALL MATMULT 
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3. Subroutine Subprogram 


SUBROUTINE ISHTAR (Y,Z) 
COMMON/1/X(100) 

Z = 0. 

DO 5 1=1,100 
5 Z = Z+X(I) 

CALL Y 
RETURN 
END 

Calling Program Reference 

COMMON/l/A(100) 
EXTERNAL PRNTIT 


CALL ISHTAR (PRNTIT,SUM) 


5.7 

A program may be written without references to subprograms or functions. On the 
other hand, it may refer to subroutines or functions or both. If so, the program is 
known as the main program. In either instance the first statement must be of the 
form: 

PROGRAM name 

where name is an alphanumeric identifier. 

A main program may refer to a variety of subroutines and functions. The sub¬ 
programs so referred to may be compiled with the main program or independently 
of the main program. Subprograms compiled with the program are called internal 
subprograms and give rise to the new terms global and local. 

An internal subprogram may be compiled with a main program, a subroutine sub¬ 
program, or a function subprogram. The first statement in the respective cases 
must be PROGRAM, SUBROUTINE S or SUBROUTINE S (p a , . . .pj, or FUNCTION 

F (pj..p n ), where the elements of the statement have established definitions. 

The method of reference to function or subroutine subprograms is the same as 
stated earlier in this chapter.. 

FORTRAN-63 assumes that all statements appearing between a PROGRAM, 
SUBROUTINE or FUNCTION statement and an END statement belong to one program. 
A typical arrangement of a set of programs and subprograms follows. 

PROGRAM SOMTHING 


END 



MAIN PROGRAM 
AND SUBPROGRAMS 
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SUBROUTINE SI 


'' END 

^ SUBROUTINE S2 
^ END 

^ FUNCTION FI (. . .) 

" END 

^ FUNCTION F2 (. . .) 

^ END 

Identifiers that are available to the entire set of programs are called global 
identifiers. Identifiers that are available only to a particular internal subprogram 
are called local identifiers. 


5.7.1 

ENTRY STATEMENT This statement provides alternate entry points to a function or subroutine 

subprogram. Its form is 

ENTRY name 

where name is an alphanumeric identifier, and may appear within the subprogram 
only in the ENTRY statement. Up to 19 entries are permitted to a given subprogram 
by use of this statement, but each entry identifier must appear in a separate ENTRY 
statement. The formal parameters, if any, appearing with the FUNCTION or 
SUBROUTINE statement do not appear with the ENTRY statement. The ENTRY 
statement may appear anywhere within the subprogram. 

In the calling program, the reference to the entry name is made just as if reference 
was being made to the FUNCTION or SUBROUTINE in which the ENTRY is imbedded. 
Rules FS4 and FS4A of 5.3.1 apply. 


Examples 

FUNCTION JOE(X,Y) 
10 JOE = X+Y 
RETURN 
ENTRY SAM 
IF (X .GR. Y) 10,20 
20 JOE = X-Y 
END 
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This could be called from the main program as follows: 


Z = A+B-JOE (3.*P,Q-1.) 


R=S+SAM(Q,2.*P) 


5.7.2 

NON-RECURSIVENESS 

OF SUBPROGRAMS Subprograms may be called from a main program or from other subprograms. Any 

subprogram called, however, may not call the calling program. That is, if program 
A calls subprogram B, subprogram B may not call program A. Furthermore, a 
program or subprogram may not call itself. 

5.8 

VARIABLE 
DIMENSIONS AND 

SUBPROGRAMS In many subprograms, especially those performing matrix manipulation, the 

programmer may wish to vary the dimension of the arrays each time the sub¬ 
program is called. This is accomplished by specifying the array identifier and 
its dimensions as formal parameters in the function or subroutine statement heading 
a function or subroutine subprogram. In die subroutine call from the calling 
program, the corresponding actual parameters are specified, and these values are 
used by the called subprogram. 


5.8.1 

The formal parameters representing the array dimensions must be simple integer 
variables in a DIMENSION statement within the subprogram. The array identifier 
must also be a formal parameter. 

VAR2 The actual parameters representing the array dimensions may be integer constants, 
integer variables, or integer arithmetic expressions. 

VAR3 If the total number of elements of a given array in the calling program is N, then 
the total number of elements of the corresponding array in the subprogram should 
not exceed N. 


VARIABLE DIMENSION 
RULES VAR1 


Example 

1. Consider a simple matrix add routine written as a subroutine subprogram: 


SUBROUTINE MATADD (X,Y,Z,M,N) 
DIMENSION X (M,N), Y(M,N), Z(M,N) 
DO 10 1=1,M 
DO 10 J=1,N 
10 Z(I,J) = X(I,J) + Y(I,J) 

RETURN 

END 
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The arrays X, Y, Z and the variable dimensions M, N must all appear as formal 
parameters in the SUBROUTINE statement and also appear in the DIMENSION 
statement as shown. If the calling program contains the array allocation 
declaration: 


DIMENSION A(10, 10), B(10,10), C(10,10) 


the program may call the subroutine MATADD from several places within the main 
program, varying the array dimension within MATADD each time as follows. 


CALL MATADD 


CALL MATADD 


CALL MATADD 


CALL MATADD 


(A, B, C, 3, 6) 


(A, B, C, 4, 8) 


(A, B, C, 3,12) 


(A, B, C, 4 * LIM, LIM2 + 3) 


When the actual parameters representing the array dimensions are integer ex¬ 
pressions, the limits of the array established by the DIMENSION statement in the 
main program may be exceeded. This condition is not checked by the compiler. 


2. Consider the 4 by n matrix: 


yn ■ 

• ■ yi.n 

r-i 

II 

• • y 2 .n 

y 3i • 

■ • y 3n 

y 4i ■ 

• • y 4n 

Its transpose Y 1 is: 


y li 

y 2 i y 3i 


Y 1 = 


^ln y2n y 3n y4n 

The following FORTRAN-63 program permits variation of n from call to call: 


SUBROUTINE MATRAN (Y, YPRIME, N) 
DIMENSION Y (4, N), YPRIME (N, 4) 

DO 7 1=1, 4 
DO 7 J= 1, N 

7 YPRIME (I, J) = Y (J, I) 

END 








% 
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APPENDIX A 


CODING FORM 


STATEMENTS 


STATEMENT SEPARATOR $ 


COMMENT CARD 


STATEMENT IDENTIFIERS 


CODING PROCEDURES 


FORTRAN-63 forms contain 80 columns in which the characters of the language 
are written, one character per column. 


The statements of FORTRAN-63 are written in columns 7 through 72. Statements 
longer than 66 columns may be carried to the next line by using a continuation 
designator. Statements may be compacted, several to a given line. Blanks may 
be used freely in FORTRAN statements to provide readability. Blanks are 
significant, however, in H fields. 



The special character $ may be used to write more than one statement on a line. 
Statements so written may also use the CONTINUATION feature. 

These statements are equivalent: 

I = 10 $ JLIM = 1 $ K = K+l $ GO TO 10 


DO 1 1=1, 10 $ A(I)=B(I)+C(I) 
1 CONTINUE $ 1=3 

1=3 


Comment information is designated by a C in column 1 of a statement. Comment 
information will appear in the source program, but it is not translated into 
object code. Columns 2 through 80 may be used. Continuation is not permitted; 
that is, each line of comments must be preceded by the column 1 C designator. 

All comment cards belonging to a specific program, or subprogram, should 
appear between the PROGRAM, SUBROUTINE, or FUNCTION statement and the 
END statement. 


I = 10 
JLIM = 1 
K = K+l 
GO TO 10 

Also: 

DO 1 1=1, 10 
A(I)=B(I)+C(I) 
1 CONTINUE 


Any statement may have an identifier (tag, label, number) but only statements 
referred to elsewhere in the program require identifiers. A statement identifier 
is a string of from 1 to 5 digits occupying any column positions 1 through 5. 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 50 59 CO 61 62 63 64 65 66 67 66 69 70 71 72 73 74 75 76 77 76 79 80 
T, A . n f£? 5081 


Hollerith Card 


1604 FORTRAN CODING FORM 


PROGRAM 


ROUTINE 




FORTRAN STATEMENT 

O = ZERO 1 

= ONE 

2 = TWO 

0 = ALPHA 0 I 

= ALPHA I 

i -- ALPHA Z 


SERIAL 

NUMBER 



FORTRAN Coding Form 



























CONTINUATION 


IDENTIFICATION FIELD 


PUNCHED CARDS 


MAGNETIC TAPE 


If I is such an identifier, 1- I - 99999; lead zeros are ignored, 1^01=001 = 0001. 
Zero is not a statement identifier. In any given program or subprogram each 
statement identifier must be unique. If the statement identifier is followed by a 
character other than zero in column 6 the statement identifier is ignored. 


The first line of every statement must have a blank in column 6. If statements 
occupy more than one line, all subsequent lines must have a FORTRAN character 
other than blank or zero in column 6. A FORTRAN-63 statement may have up to 
600 alphanumeric characters, operators, delimiters (commas or parentheses) 
and identifiers within it; blanks are not included in this count. In general, up to 
20 continuations may appear after a first statement. 


Columns 73 through 80 are always ignored in the translation process. These 
columns, therefore, may be used for identification when the program is to be 
punched on cards. Usually these columns contain sequencing information 
provided by the programmer. 


Each line of the coding form corresponds to one 80-column card, and the terms 
line and card are often used interchangeably. Source programs and data can be 
read into the computer from cards; an object program memory map, or data, 
can be punched directly onto cards. Usually, however, cards are used in the 
off-line preparation of input or output magnetic tapes. 

Blank cards appearing within the input card deck are treated as follows: 

a) If a blank card appears between a statement and its continuation, the 
continuation and other continuations following it are lost. Compilation 
continues. 

b) If a blank card appears between two statements, it is ignored. 

When cards are being used as a data medium rather than as statements, all 80 
columns may be used. 


Magnetic tapes are the most commonly used input-output media. Tape charac¬ 
teristics are described in the Control Data Corporation hardware manuals. The 
record structure resulting from the READ/WRITE or BUFFER I/O control 
statements is described in Chapter 2 of Volume II. 

Compilation of FORTRAN-63 programs requires the following tapes: 

Master tape 
Standard mput unit 
Standard output unit 

FORTRAN scratch tape I not used for very 
Assembler scratch tape / small programs 
Standard punch unit 
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CARRIAGE CONTROL 


The first character of the printer record is a control character for providing line 
spacing control. This character is not printed but will cause the following actions: 


Character 

Blank 

0 

1 


Action 

Single space before printing 
Double space before printing 
Eject page before printing 


These codes are standard on all printers used with the 1604. Some printers 
provide additional codes which are given in the specific manuals. 



APPENDIX B 


CHARACTER CODES 1604 COMPUTER 


Source Language 
Character 


BCD (Magnetic Punch Positions in a 

Tape & Internal) Hollerith Card Column 


A 

61 


12-1 

B 

62 


12-2 

C 

63 


12-3 

D 

64 


12-4 

E 

65 


12-5 

F 

66 


12-6 

G 

67 


12-7 

H 

70 


12-8 

I 

71 


12-9 

J 

41 


11-1 

K 

42 


11-2 

L 

43 


11-3 

M 

44 


11-4 

N 

45 


11-5 

0 

46 


11-6 

P 

47 


11-7 

Q 

50 


11-8 

R 

51 


11-9 

S 

22 


0-2 

T 

23 


0-3 

U 

24 


0-4 

V 

25 


0-5 

W 

26 


0-6 

X 

27 


0-7 

Y 

30 


0-8 

Z 

31 


0-9 

0 

12 


0 

1 

01 


1 

2 

02 


2 

3 

03 


3 

4 

04 


4 

5 

05 


5 

6 

06 


6 

7 

07 


7 

8 

10 


8 

9 

11 


9 

/ 

21 


0-1 

+ 

60 


12 

- 

40,14 


11,8-4 

lank 

20 


space 


73 


12-8-3 

) 

74 


12-8-4 

$ 

53 


11-8-3 

* 

54 


11-8-4 


33 


0-8-3 

( 

34 


0-8-4 

= 

13 


8-3 
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CHARACTER CODES 3600 COMPUTER 


Source Language 

BCD (Internal 

Punch position in a 

Character 

only) * 

Hollerith Card Column 

A 

21 

12-1 

B 

22 

12-2 

C 

23 

12-3 

D 

24 

12-4 

E 

25 

12-5 

F 

26 

12-6 

G 

27 

12-7 

H 

30 

12-8 

I 

31 

12-9 

J 

41 

11-1 

K 

42 

11-2 

L 

43 

11-3 

M 

44 

11-4 

N 

45 

11-5 

O 

46 

11-6 

P 

47 

11-7 

Q 

50 

11-8 

R 

51 

11-9 

s 

62 

0-2 

T 

63 

0-3 

U 

64 

0-4 

V 

65 

0-5 

w 

66 

0-6 

X 

67 

0-7 

Y 

70 

0-8 

Z 

71 

0-9 

0 

00 

0 

1 

01 

1 

2 

02 

2 

3 

03 

3 

4 

04 

4 

5 

05 

5 

6 

06 

6 

7 

07 

7 

8 

10 

8 

9 

11 

9 

/ 

61 

0-1 

+ 

20 

12 

- 

40 

11-8-4 

blank 

60 

space 


33 

12-8-3 

) 

34 

12-8-4 

$ 

53 

11-8-3 

* 

54 

11-8-4 

, 

73 

0-8-3 

( 

74 

0-8-4 

= 

13 

8-3 


^Magnetic Tape Codes same as 1604. 
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APPENDIX C 


REPLACEMENT 


CONTROL 


TYPE 

DECLARATION 


* E = Executable 


STATEMENTS OF FORTRAN-63 


Page 


A=E Arithmetic 

E* 

14 

A=L Logical 

E 

20 

A=M Masking 

E 

21 

.=Aji=E Multiple 

E 

22 

GO TO n 

E 

37 

GO TO n,(n!, . . .n m ) 

E 

37 

GO TO (n x , . . .,n m ),i 

E 

38 

ASSIGN i to n 

E 

37 

A)n! ,n 2 ,n 3 

E 

38 

IF( b )nx ,n 2 

E 

39 

SENSE LIGHT i 

E 

39 

IF(SENSE LIGHT i)m ,n 2 

E 

39 

IF(SENSE SWITCH i)m ,n a 

E 

39 

(FAULT 1 

E 

39 

IF DIVIDE { CHECK } s n 2 

IF EXPONENT FAULT n x ,n a 

E 

39 

IF OVERFLOW FAULT n x ,n a 

E 

40 

DO n i=m x ,m 2 ,m 3 

E 

40 

CONTINUE 

E 

42 

PAUSE 

E 

43 

STOP 

E 

43 

END 

N/E 

43 

TYPE COMPLEX List 

N 

25 

TYPE DOUBLE List 

N 

25 

TYPE REAL List 

N 

25 

TYPE INTEGER List 

N 

25 

TYPE LOGICAL List 

N 

25 

TYPE named (W/f) List 

N 

25 


d is 5, 6, or 7 


N = Non-executable 




Volume I 
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Statements of FORTRAN-63 (Continued) 


STORAGE 

DIMENSION Vj ,V 2 , ... 

N 

Page 

32 Volume I 

ALLOCATION 

COMMON/B(I)/List . . . 

N 

27 


EQUTVALENCE(a,b,c. . .)(p,q. . .) 

N 

30 

DATA 

DATA(I=List) ,(I=List), ... 

N 

32 

STATEMENT 

SUBPROGRAM 

FUNCTION name(p. ,p 2 ,_) 

N 

45 

STATEMENTS 


SUBROUTINE name(p 1 ,p 2 , . . .) 

N 

46 


PROGRAM name 

N 

47 


EXTERNAL namej ,name 2 , . . . 

N 

48 


ENTRY name 

N 

53 


CALL name 

E 

48 


RETURN ^ 

E 

48 

I/O FORMAT 

FORMAT(spec 1 ,spec 2 , . . .) 

N 

Volume n 

I/O READ/WRITE 

READ n,I 

E 

25 


PRINT n,L 

E 

21 


PUNCH n,L 

E 

21 


READ(i,n)L ) 

E 

21 


READ INPUT TAPE i,n,L f 

WRITE(i,n)L ) 

E 

21 


WRITE OUTPUT TAPE i,n,L ( 

READ(i)L \ 

E 

25 


READ TAPE i,L ) 

WRITE(i)L J 

E 

25 

I/O TAPE 

WRITE TAPE i,L ) 

END FILE i 

E 

26 

HANDLING _ 


REWIND i 

E 

26 


BACKSPACE i 

E 

26 

I/O STATUS 

IF (EOF.i)^ ,n 2 

E 

28 

CHECKING 

IF (IOCHECK.i)^ ,n 2 

E 

28 


IF (UNIT,i)n 1 ,n 2 ,n 3 ,n 4 

E 

28 

I/O BUFFERING 

BUFFER IN(i,p)(A,B) ) 

E 

26 

INTERNAL DATA 

BUFFER OUT (i,p)(A,B) ( 

ENCODE (c,n,V)L 

E 

29 

MANIPUATION __ 


DECODE (c,n,V)L 

L 

2y 



APPENDIX D 


LIBRARY FUNCTIONS 


Form 

Definition 

Actual Para¬ 

Mode of 



meter type 

Result 

ABSF(X) ( 


Real 

Real 

XABSF (i)INTF (X) ( 

Absolute Value 

Integer 

Integer 

INTF(X) \ 


Real 

Real 

XINTF(X) ( 

Truncation, integer 

Real 

Integer 

MODF(X 1 ,X 2 ) 

X, modulo X 2 

Real 

Integer 

XMODF^d i 2 ) 

ij modulo i 2 

Integer 

Integer 

MAX0F(i l9 i 2 . . . ) 

\ 

Integer 

Real 

MAX1F(X 1} X 2 , . . . ) 

/ 

Real 

Real 

XMAX0F(i l9 i 2 , . . . ) 

>Determine maximum argument 

Integer 

Integer 

SMAX1F(X 15 X 2 , . . . 

)) 

Real 

Integer 

MIN0F(i 1} L,, . . . ) 

\ 

Integer 

Real 

MINlFCXj,X 2 , . . . ) 

1 

Real 

Real 

XMIN0F(i l5 i 2 , . . • ) 

/Determine minimum argument 

Integer 

Integer 

xminif(x 15 x 2 , . . . : 

>) 

Real 

Integer 

SINF (X) 

Sine X radians 

Real 

Real 

COSF(X) 

Cosine X radians 

Real 

Real 

TANF(X) 

Tangent X radians 

Real 

Real 

ASINF(X) 

Arcsine X radians 

Real 

Real 

ACOSF(X) 

Arccos X radians 

Real 

Real 

ATANF(X) 

Arctangent X radians 

Real 

Real 

TANHF(X) 

Hyperbolic tangent X radians 

Real 

Real 

SQRTF(X) 

Square root of X 

Real 

Real 

LOGF(X) 

Natural log of X 

Real 

Real 

EXPF(X) 

e to xth power 

Real 

Real 

SIGNF(X 1s X 2 ) 

Sign of X 2 times X, 

Real 

Real 

XSIGNF(ii ,i 2 ) 

Sign of i 2 times ij 

Integer 

Integer 

DIMF(X x ,X 2 ) 

( If Xj > X 2 : Xj-Xa 




\ 

Real 

Real 


{ If XI ^ X 2 : 0 



XDIMF(ii ,1.2) 

{ If i x > i 2 : ii -i 2 




(If i x ^ i 2 : 0 

Integer 

Integer 
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Form 

Definition 

Actual Para¬ 

Mode of 



meter Type 

Result 

CUBERTF (X) 

Cube root of X 

Real 

Real 

FLOATF(I) 

Integer of Real Conversion 

Integer 

Real 

RANF(N) 

Generate Random Number 

-Real ( 

-Integer j 

$ 

Real 


(Repeated Executions give 

+Real | 

Integer 


uniformly distributed numbers) 

+Integer ) 

XFIXF 

Real to Integer Conversion 

Real 

Integer 

POWERF(X 15 X 2 ) 

M X 

to 

Real, Real 

Real 

ITO J (I, J) 


Integer, Integer 

Integer 

XTOI(X,I) 

X 1 

Real, Integer 

Real 

ITOX(I,X) 

I s 

Integer, Real 

Real 
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INDEX - Volume I 


Arithmetic Expressions 

Page 

9 

Arithmetic Replacement Statement 

14 

Array Structure and Subscripts 

6 

ASSIGN Statements 

38 

Assigned GO TO 

37 

Characters 

1 

Codes 

62 

1604 Character 

62 

3600 Character 

63 

Coding Procedures 

58 

$ 

Comments 

58 

COMMON 

27 

COMMON Blocks 

28 

Computed GO TO 

38 

Constants 

3 

Hollerith 

3 

Integer 

3 

Octal 

3 

CONTINUE 

42 

Control Statements 

37 

Conversions, Mixed Mode 

14 

DATA 

32 

DIMENSION 

2,6 

Divide Check/Fault 

39 

DO Statement 

40 

END 

43, • 

ENTRY 

53 

EQUIVALENCE 

30 

Exponent Fault, If 

39 

Expressions 

9 

Arithmetic 

9 

Logical 

16 

Masking 

20 

Mode of Arithmetic 

12 


EXTERNAL 
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INDEX - Volume I 

Page 


Fault Condition Statements 39 

If Divide Check/ 39 

If Exponent 39 

If Overflow 40 

FORTRAN-6 3 Language Statements 8 

Function Subprograms 46 

Functions 45 

Library 46,66 

Statement 45 

GO TO Statements 37 

Assigned 37 

Computed 38 

Unconditional 37 

Hollerith Constants 3 

Identifiers 2 

Statement 37,58 

IF Statements 38 

Divide Check/Fault 39 

Exponent Fault 39 

Overflow Fault 40 

(SENSE LIGHT i) 39 

(SENSE SWITCH i) 39 

Three Branch 38 

Two Branch 39 

Integer Constants 3 

Library Functions 46 

Logical Expressions 16 

Logical Replacement Statement 20 

Main Program and Subprograms 52 

Masking Expressions 20 

Masking Replacement Statement 21 

Mode of Arithmetic Expressions 12 

Multiple Replacement Statement 22 
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INDEX - Volume I 




Octal Constants 

Page 

3 

Operators 

1 

Overflow Fault 

40 

PAUSE 

43 

Quantities Structure 

2 

Replacement Statement 


Arithmetic 

14 

Logical 

20 

Masking 

21 

Multiple 

22 

RETURN 

48 

Statement Functions 

45 

Statement Identifiers 

37, 58 

Statement Separator 

58 

Statements of FORTRAN-63 

64 

STOP 

43 

Storage Allocation 

25 

Subroutines 

45 

Subprograms 

50 

Type Declaration 

25 

Variables 

6 

Dimensions 

54 

Subprograms 

54 

Word Structure 

2, 5 


i 
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Errata Sheet 


FORTRAN-63/Reference Manual - Volume I continued - Publication Number 527 


>age 

Article 

Remarks 

59 

Appendix A 

FORTRAN Coding Form, Statement 30 should read 
"LOGR = LOGF (R) 11 

63 

Appendix B 

Punch position for minus is 11,8-4 

66 

Appendix D 

SMAX1F should be XMAX1F 
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Errata Sheet 



FORTRAN- 

-63/Reference Manual Volume 1 Publication #527 

Page 

Article 


17 

2.5 

First line: strike the character "/". 

27 

3.3 

Last paragraph before examples. Add the sentence "Leading 
zeros in numeric identifiers are ignored." 

32 

3.5.2 

Second paragraph; change "If any, or all," to "If all". 

33 

3.6.1 

Rule DAT6, 2nd line should read in part: 

"...DATA (A=2), the type ...” 

39 

4.3.3 

The F63 statement SENSE LIGHT is omitted. It reads: 

"SENSE LIGHT i 

The statement turns on the i-th sense light 0 SENSE LIGHT 0 
turns off all sense lights. i may be a simple integer 
variable or constant. 

In the 3600, 1 £, i ^ 48 

In the 1604, 1 £t i ^ 4<" 

39 

4 0 3.4 

Strike the words "Appendix E". 

40 

4.5 

Strike out the last sentence of the last paragraph. Strike 
out the sentence "The DO statement provides FORTRAN-63, etc. 

46 

5.1 

Rule SF3. Add: "Formal parameters must be simple variables 

Rule SF5. Change to: "E may contain subscripted variables, 
but the subscripts are restricted to integer constants." 

47 

5.3.1 

Rule FS4, next to last sentence--last word should be "list", 
not "bit". 

49 

5.5 

Example 4: 

"END" follows "RETURN" 

53 

5.7.1 

In the Example: 

Change "ENTRY SAM" to "ENTRY JAM" 

Change "IF (X.GR.Y)10,20" to "IF(X.GT.Y)10,20". 

54 

5.7.1 

Change "S + SAM(Q,2.*p)" to "S + JAM(Q,2.*P)" 

54 

5.8.1 

Add to rule VAR1: 

"The formal parameters must not appear in a COMMON or 
EQUIVALENCE statement in the subprogram”. 

55 

5.8.1 

Example 2 


Change "DO 7 1=1,4 to "DO 7 1=1,N 
DO 7 J=1,N" DO 7 J=l,4" 
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